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We show that violation of the null energy condition implies instability in a broad class of models, 
including classical gauge theories with scalar and fermionic matter as well as any perfect fluid. When 
applied to the dark energy, our results imply that w = p/p is unlikely to be less than — 1. 



I. INTRODUCTION 

Energy conditions, or restrictions on the matter 
energy-momentum tensor T^ v , play an important role 
in general relativity. No classification of the solutions 
to Einstein's equation is possible without restrictions on 
T M „ , since every spacetime is a solution for some particu- 
lar choice of energy-momentum tensor. In this letter we 
demonstrate a direct connection between stability and 
the null energy condition (NEC) Q, T^n^n" > for 
any null vector n (satisfying g^n^n 1 ' — 0). Our main 
results are: (1) classical solutions of scalar-gauge mod- 
els which violate the NEC are unstable, (2) a quantum 
state (including fermions) in which the expectation of 
the energy-momentum tensor violates the NEC cannot 
be the ground state, (3) perfect fluids which violate the 
NEC are unstable. These results suggest that violations 
of the NEC in physically interesting cases are likely to be 
only ephemeral. 

Our results have immediate applications to the dark 
energy equation of state, often given in terms of w = pj p. 
Dark energy has positive energy density p and energy- 
momentum tensor T M „ — diag (p, p, p, p) in the comoving 
cosmological frame. Therefore, w < — 1 implies violation 
of the NEC. Instability as a consequence of w < —1 was 
studied previously in scalar models Q. 

Some results in relativity in which the NEC plays 
an important role include the classical singularity theo- 
rems |3j , proposed covariant entropybounds |J] and non- 
existence of Lorcntzian wormholcs 5]. 



II. FIELD THEORIES 

Consider a theory of scalar, <j) a , and gauge, A aa , fields 
in a fixed c?-dimensional space-time with the metric g^ v . 
We limit ourselves to theories whose equations of motion 
are second order differential equations, so the Lagrangian 
for the system is assumed to depend only on the fields 
and their first derivatives. We take the Lagrangian den- 
sity C to depend only on the covariant derivative of the 
field D^ipa and the gauge field strength F afiv . The scalars 
may transform in any representation of the gauge group. 
We impose Lorentz invariance on C, but do not require 



overall gauge invariance. That is, we allow for fixed ten- 
sors with gauge indices (but no Lorentz indices) which 
can be contracted with the fields. For the corresponding 
action 



S = 



J d d x\g\i£{(j) a ,D^ a ,F a ^) 



(1) 



to be stationary, its first variation has to vanish, SS = 0. 
This leads to the equations of motion for the fields <j) a 
and A aa ; in the classical analysis we assume that we 
have found solutions to these equations, about which we 
expand. 



A. Null energy condition 

The quantities D^<j) a and g^ v are independent vari- 
ables. Nevertheless we now prove that there is a relation 
between the derivatives of C with respect to them: 



= M AB ip^ 



(2) 
(3) 
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In our notation ipA^ = {Dn4>a, F aa ^), where the abstract 
index A may run over both Lorentz and color indices, as 
well as the type of field. So, tpAn is a list of objects, each 
of which has a Lorentz index fi. The value of A specifies 
an element of this list. 

The relations are obtained by noting that for each 
and every g^ v in C there are two ^s attached to it, ex- 
cept for the curved space totally antisymmetric tensor 
\g\ ~ Sf"! •••"d i which gives rise to the K term in Eq. J2J. 
Similarly, differentiation with respect to ipAn yields the 
M and L terms in Eq. ©. 

Figure ^ represents the most general Lagrangian of 
type considered in this paper. Each dot represents a 
Lorentz index and a line connecting them denotes con- 
traction using the metric. Rectangles (with two indices) 
are field strengths, small circles covariant derivatives of 
scalar fields, and a large circle an epsilon tensor. Finally, 
the block X represents the remainder of the diagram. Be- 
cause the product of two epsilon tensors can be rewritten 
as a sum of products of metric tensors g, we need to con- 
sider only figures with at most one epsilon tensor, and 
therefore can assume that X contains none. (For gener- 
ality, we include an epsilon tensor in the figure, although 
of course C need not contain a parity- violating compo- 
nent.) All Lorentz indices are ultimately contracted, and 
we suppress color indices for simplicity. 
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FIG. 1: Representation of the most general Lagrangian of 
the type considered in this paper. Each dot represents a 
Lorentz index, and a line connecting them denotes contrac- 
tion using the metric. Rectangles (with two indices) are field 
strengths, small circles covariant derivatives of scalar fields, 
and large circles epsilon tensors. Finally, the block X repre- 
sents the remainder of the diagram. All Lorentz indices are 
ultimately contracted, and we suppress color indices for sim- 
plicity. In graphical terms, M is obtained by simply removing 
the shaded elements. 



Consider the labelled portion of the figure, which 
equals 



1.91 



(4) 



Again, we suppress color indices for simplicity, as they do 
not affect the proof. We include an epsilon tensor in the 
analysis, although £ may or may not contain one (in the 
parity-preserving case the epsilon tensor in Eq. (0J is re- 
placed by the metric) . By differentiation of the indicated 
portion, we obtain Eqs. @ and with 



M af3 
K 
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(5) 
(G) 
(7) 



Note the derivative £f„ i3 generates a contribution to M 
which is matched by a corresponding contribution from 
2£ 9 /3 T . Other contractions of fields with g pl/ (i.e., as indi- 
cated in the figure) can be analyzed similarly. In graph- 
ical terms, M can be obtained from the figure for £ by 
simply removing two ?/>s, in this case the shaded elements 
of the figure. 

For the energy-momentum tensor which couples to 
gravity, 



T^ u £g -\- 2£g^^ , 
the NEC then requires 

^ A M AB ^ B > 0, 



(8) 



(9) 



rAB 



where * A = ipA^. Thus, to satisfy the NEC, AP 
has to be positive semidefinite. This property is crucial 
for stability of solutions, to which we now turn. 



B. Stability 

To study the stability of the solution ipA (x) , we con- 
sider the second variation of the Lagrangian, 

6 2 £ = C^ A1 p B 5ijjAHB + 2C^ a ^ b . x SiIjaS'4>b;X 

+ ^a-^-b^Ha-^Hb-v (10) 

Here quantities £^ A = d£/dipA, etc. are evaluated at 
■>Pa{x). Also notice that ipA-^i = (D^a, A oce;/1 ), the 
covariant derivatives of ipA, are different from ipAn — 
(FJ^cfya, F aap ). 

Let us use a locally inertial frame in which the met- 
ric is reduced to g^ = diag (1, — 1, . . . , — 1); all quan- 
tities in this frame are designated by a bar. For the 
Lagrangian i|l(J|) . the canonical momentum is 

8H A = 2£„ 



Mi 



2£ 



i>A;0"<Pl 



Sip 



B;u 



(11) 



which leads to the following effective Hamiltonian for 
fluctuations about the classical solution: 

8 2 H = -L^ a ^J^aHb-2L^ a ^ b J^aHb^ (12) 

Here 6tjjA;0 are functions of Stt b , 5ipB and 5ipB;i as found 
from Eq. (|11|) . The first term on the right hand side of 
Eq. (|12fl is a potential term, which we denote by 5 2 V, 
and the last two terms are kinetic terms, denoted S 2 IC. 

If the kinetic energy 5 2 JC is negative, then the system 
described by the Hamiltonian 8 2 TL = 8 2 JC + 8 2 V is (lo- 
cally) unstable. If 8 2 V is positive then small perturba- 
tions will cause the classical solutions to grow exponen- 
tially away from the original stationary point. However, 
it is possible to have classical stability if one chooses 8 2 V 
to be negative; in this case we have an upside-down po- 
tential with negative kinetic term, or a "phantom" . Such 
models necessarily exhibit quantum instabilities No- 
tice, the second term in Eq. I|12|) is linear in the fluc- 
tuations and their derivatives, and therefore can never 
stabilize the system. 

To investigate the kinetic terms, we calculate second 
derivatives of £ needed in Eq. (fT^fl . Using Eq. 10) we 
obtain 



C MB „=M AB g^ + N A » 



Br 



(13) 



The separation of the second derivative into the first and 
second terms in Eq. is natural: g^ appears only 
in the first term, and N represents all remaining terms. 
A" contains terms obtained by differentiating M and L 
with respect to ipBv, plus additional terms if tp is a field 
strength. The v index obtained from these ipBu deriva- 
tives is attached to a field and not the metric g^ . (Also, 
L does not contain an epsilon tensor since X does not.) 
Finally, notice that even though ipAp, and ipA-^ differ, the 
derivatives of £ with respect to them coincide due to the 
form of the action J[J . Thus the kinetic term becomes 

(M AB + N AQB0 ) 5j A flS$ B .o 
(M AB S*3 - N MB >) StPa-Mb-j. (14) 



S 2 JC 



3 



We now prove that nonnegativeness of the kinetic term if 



6 JC implies positive semidefiniteness of the matrix M. 
Indeed, suppose that M is not positive semidefinite. 
In such case, the matrix M has at least one negative 
eigenvalue, which means that there is a basis in which 
the matrix is diagonal with at least one negative entry, 
M = diag (mi, . . . , m n ) (mi < 0). (Quantities in this ba- 
sis are designated with a tilde.) Let us choose such field 
variations that are nonzero only in the direction of the 
negative eigenvalue: S-ipi^ ^ and StpA-^ = (A > 1). 
We further restrict d — 1 quantities dipi-i to satisfy the 
following equation: 



(15) 



These conditions make the kinetic term of Eq. (|14|) neg- 
ative, 



S JC = mi 



<0, 



(16) 



thus proving that in order for 5 2 JC to be nonnegative, the 
matrix M has to be positive semidefinite. 

Using the result established in the previous paragraph, 
we conclude that solutions to the theory given by the 
action (JIJ are stable only if the matrix M is positive 
semidefinite. 

Combining the relations between nonnegativeness of 
S 2 IC and positive semidefiniteness of M on one hand, and 
the NEC and positive semidefiniteness of M on the the 
other hand, we conclude that for the theory given by 
the action (JIJ, only solutions satisfying the NEC can be 
stable. 

We can deduce similar results for quantum systems. 
Suppose there exists a quantum state \a) and a null vec- 
tor n M such that 



{a\T liV \a)ni>n v = {a\M A * V A V B \a) < 



(17) 



so that the NEC is violated in a quantum averaged 
sense. Define a basis \(j>) in which the operator M = 
M ab ^ a ^b is diagonal: 

M\^)=m{4>)\4>). (18) 

Then violation of the quantum averaged NEC implies 

Y^WmMWWW) = Y, K«l$l 2 mfa)<0. (19) 

</></>' <t> 

This means that there exist eigenstates \<j>), whose over- 
lap with \a) is non-zero and on which the operator M. 
has negative eigenvalues. This requires that M and hence 
S 2 IC is not positive semidefinite; by continuity, this must 
also be the case in a ball B in the Hilbert space of \<f>). 

As a further consequence, we can conclude that a 
state | a) in which the NEC is violated cannot be the 
ground state @- Suppose that \a) is an energy eigen- 
state: H\a) — E a \a). An elementary result from quan- 
tum mechanics is that \a) can be the ground state only 



E a = (a\H\a) < {a'\H\a') 



(20) 



for all normalized states \a') which need not be energy 
eigenstates. However, it is possible to reduce the expecta- 
tion value of H by perturbing \a). Specifically, we adjust 
\a) only in the ball B, where we know from Eqs. (|14|) — (|lfc>|> 
that there are perturbations which reduce the expecta- 
tion of the kinetic energy without changing the expecta- 
tion of the potential. 

Note that the discussion above is in terms of unrenor- 
malized (bare) quantities. The renormalized expectation 
(a\M Ten \a) = (a\M\a) - (0\M\0) (where |0) is the flat- 
space QFT ground state) could be negative (e.g., as in 
the Casimir effect pj), but this is possible only if \a) is 
not |0). 

In known cases of NEC violation, such as the Casimir 
vacuum or black hole spacetime, it is only the renormal- 
ized energy-momentum tensor which violates the NEC. 
As a simple example, consider a real scalar field cf>. The 
energy-momentum tensor is simply T M „ — d^ipd^ip plus 
terms proportional to g^ v which do not play a role in 
the NEC. Then, M — (n^d^cj)) 2 is a Hermitian opera- 
tor with positive eigenvalues. Therefore, its expectation 
value in any state is positive: (a|.A/(|a} > 0, for any |a), 
including the Hartlc-Hawking, Casimir or flat-space vac- 
uum. We can verify this by direct calculation, computing 
the energy-momentum tensor using point-splitting regu- 
larization: 



(21) 



which is manifestly positive. Note that (a\A4\a) > for 
all | a), since the bare expectation is always dominated 
by the UV contribution. Now, had we taken a negative 
kinetic energy term for the scalar, the overall sign of Ai 
would change, allowing violation of the NEC. But, this 
model is clearly unstable, in accordance with our results. 



C. Fermions 

To this point we have only considered bosonic fields. 
We now extend our analysis to systems with fermions, 
adding to our Lagrangian the term 



(22) 



(A scalar-fermion coupling can be treated similarly, as 
can a Weyl fermion, whose determinant is the square root 
of the Dirac determinant.) Then, for any fixed gauge field 
background the fermions can be integrated out directly 
in favor of a non-local correction to the action for bosonic 
fields 



E 

Ai>0 



ln(A 



(23) 
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where ftipi = Xiipi is the eigenvalue equation for the Dirac 
operator, with A; real. This shifts the energy-momentum 
tensor by 

T$ = -2\g\-i Y N9 1 ? f ' . (24) 

Now write A 2 ^ = = {g^D„D v - F^)^ and 

use the orthonormality of the eigenfunctions to obtain 

X 2 = J 'd d x\g\^ g^^D^t + (25) 

where the ellipsis denote terms which do not contain g^ v . 
After integration by parts, the contribution to the NEC 
from fermions is then 

r$n"n" - ]T TTT^I (» ' ^) + (" • Dft). (26) 

A/ — r~ Til 
A,>0 < 

This additional contribution is always positive. So, the 
conclusions of the previous section are unmodified by the 
presence of fermions: violation of the NEC implies the 
bosonic kinetic energy is not positive semidefmitc. 

III. PERFECT FLUID 

A macroscopic system may be approximately described 
as a perfect fluid if the mean free path of its components 
is small compared to the length scale of interest. For the 
dark energy, this length scale is of cosmological size. A 
perfect fluid is described by the energy-momentum tensor 

Tpu = (p + p)u^u v - pg^ v , (27) 

where p and p are the energy density and pressure of 
the fluid in its rest frame, and is its velocity. Let 
j 1 * = Ju^ be the conserved current vector {j^.^ = 0), 

and J = (j^j^)? the particle density. 

The energy-momentum can be written 0, Q as 

T MV = (f- Jf')9»v + (f'/J^ju, (28) 
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